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>► ! Abstract 

We summarize results on the asymptotics of the two-particle Green 
f-x ■ functions of interacting electrons in one dimension. Below a critical 

value of the chemical potential the Fermi surface vanishes, and the 
system can no longer be described as a Luttinger liquid. Instead, the 
non-relativistic Fermi gas with infinite point-like repulsion becomes 
the universal model for the long-wavelength, low temperature physics 
of the one-dimensional electrons. This model, which we call the im- 
penetrable electron gas, allows for a rigorous mathematical treatment. 



O In particular, a so-called determinant representation for the two par- 

; ticle Green function could be derived. This determinant represen- 

tation is related to an integrable classical evolution equation and to 
a Riemann-Hilbert problem, that enable the exact calculation of the 
asymptotics of the two-particle Green functions. 
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1 Correlation functions 

Correlation functions are (thermal) expectation values of products of field oper- 
ators. As a typical example we may consider the two-particle Green function 
G a p(x,t) = (*P a (^)0^p (OjO))- Here and are annihilation and creation 
operators of electrons with spin a =f, |. The brackets denote the grand canonical 
ensemble average. The two-particle Green function G a p(jc, t) gives the probability 
to find an electron at (x, t) provided there was an electron at (0, 0) . It describes the 
propagation of an electron in a medium of other electrons, which, in general, is 
characterized by a temperature T, a chemical potential /j and a magnetic field B. 

Within the frame work of the linear response theory correlation functions pro- 
vide the link between the microscopic and the macroscopic properties of mat- 
ter. Correlation functions are measurable in experiments. The two-particle Green 
function, for instance, measures the optical absorption. 

Only a few mathematically exact results on correlation functions of interacting 
systems are known. In fact, the example discussed in this contribution is the first 
example of a direct Bethe ansatz calculation of the asymptotics of correlations of 
interacting electrons. The calculation was performed for the so-called impenetra- 
ble electron gas model, which is the infinite repulsive coupling limit of Yang's 
model of electrons with point-like pair interaction [|l|, §]. As we shall argue be- 
low, our results are, in spite of being obtained for a specific model, applicable to 
a whole class of models of interacting electrons in one dimension, and are, in this 
sense, universal. 

In the first part of this text we shall explain our ideas about the universality 
of the long-wavelength, low-temperature asymptotics of electronic correlations 
in the gas phase [JJ]. We shall start with the paradigmatic Hubbard model, and 
then argue that certain modifications of the interaction do not change the long- 
wavelength, low-temperature physics of the model. In the appropriate scaling 
limit all modified Hamiltonians lead to the same effective model. The second part 
of this text is devoted to a summary of our calculation of the asymptotics of the 
two-particle Green functions of the impenetrable electron gas model [Q, ||]. 



The density D of non-interacting, one-dimensional spin-^ Fermions on a lattice is 
given by the integral over the Fermi weight, 



Here e(p) is the dispersion of the Fermions, T denotes the temperature and fj 
the chemical potential. Let us assume z{p) to be monotonically increasing and 



2 The Hubbard model in the gas phase 




(i) 
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bounded from below. If the chemical potential is smaller than a critical value, 
p < p c = irhn p> ()£(p), then D vanishes in the zero temperature limit T — > 0+. 
For p > p c , on the other hand, the density D approaches a finite positive value as 
T — > 0+. This means that the system undergoes a phase transition at T = as a 
function of the chemical potential. The critical point is at ^ = p c . Assuming that 
e(p) = p c + p 2 + O (p 4 ), we obtain 



D=-^p7 c (2) 
n 

for p c < p < p c + 8, 8 <C 1 . 

As we shall see below by considering a representative example this scenario 
remains unchanged if the Fermions interact. For interacting one-dimensional 
Fermions the phase with p > p c , for which the density at T = is finite, is called 
the Luttinger liquid phase [|J [7|, |8|| . This phase is quite familiar to many physicists. 
A Luttinger liquid may be understood as a one-dimensional metal. The correla- 
tions in the Luttinger liquid are dominated by fluctuations around the Fermi sur- 
face. Their power law decay at zero temperature is described by conformal field 
theory . For small finite temperature conformal field theory predicts exponen- 
tially decaying correlations. One has to employ a conformal mapping from the 
complex plane to a strip of finite width. As a result the rate of exponential decay 
is defined by conformal dimensions. 

The phase with p < p c so far attracted less attention. This phase is trivial at 
T = 0, since the density D vanishes for T = 0. The density becomes positive for 
positive temperature, and is typically exponentially small as long as the tempera- 
ture remains small. The ideal gas law holds. This suggested the name 'gas phase' 
to us. The gas phase may be interpreted as a one-dimensional semi conductor or 
insulator. Correlations in the gas phase behave essentially different compared to 
those in the Luttinger liquid phase. This is the subject of this text. 

In order to get a better understanding of the gas phase of interacting systems 
let us consider the Hubbard model as an example, 

L L L 

^ = -£(c+ +La c i , CT + c+ c ;+ i, CT ) + ^£n jT n ji -A/£(« J - T + n ;i ) . (3) 

Here the canonical Fermi operators cj , Cj, a are creation and annihilation oper- 
ators of electrons at site j of a one-dimensional, periodically closed chain of L 
lattice sites, and nj j and n; i are the corresponding particle number operators. U 
is the strength of the repulsive interaction. 

The eigenvalue problem of the Hubbard Hamiltonian can be solved JTC| ] by 
means of the nested Bethe ansatz. This allows us to test our ideas about the gas 
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phase quantitatively. The energy levels for the N electron system are 



N 



E = 2 £ (1 - cos*,-) - 2)N , 



(4) 



where the charge momenta kj are solutions of the Lieb-Wu equations [{Hip . Clearly 
the first term on the right hand side of (Q) is non-negative. Hence, if /j < /j c = — 2, 
the energies of all eigenstates become non-negative, and the absolute ground state 
is the empty lattice. For /j > —2, on the other hand, the energy can be lowered by 
filling states with small fc's. Since kj+i — k; ~ 1/L, this leads to a finite density of 
electrons in the ground state as L — ► °°. We conclude that the Hubbard model is 
in the gas phase for /j < —2 and in the Luttinger liquid phase else. Note that the 
asymptotics of correlation functions of the Hubbard model in the Luttinger liquid 
phase was obtained in [ JITT , [12| ] . 

Another way of understanding the transition from the gas phase to the Lut- 
tinger liquid phase is by looking at the integral equations that describe the Hub- 
bard model in the thermodynamic limit (see e.g. JT3|]). At zero temperature and 
zero magnetic field the dressed energy K(k) of the elementary charge excitations 
is determined by the integral equation 



K(k) = -2cos(fc) -fi+ dk' cos(k')R(sm(k') - sm(k))K(k!) , (5) 



where the limits of integration depend on the chemical potential through the con- 
dition 



Similarly, the density p(k) of elementary charge excitations is obtained from the 
integral equation 




K(fi)=0. 



(6) 



The integral kernel R is given by 




(7) 




(8) 



p(&) determines the density of the electrons at zero temperature, 



D 



J-Q 



/ dkp(k) 



(9) 
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Since p(k) is positive, D is equal to zero if and only if Q = 0. Hence, it follows 
from (Q) and (§) that /j c = —2. Moreover, the equations ©-(0) are easily solved 
for small Q. From the solution we obtain 

1 



D=-^JT C (10) 
n 

for n c < /j < /j c + 8, 8 < 1 . 

Thus the qualitative picture is the same for free and for interacting electrons. 
Comparing @ and (|10|), however, we see that the results for the electron density 
close to the critical point /j = /j c differ by a factor of two. This difference may be 
interpreted as a signature of a phenomenon called spin-charge separation (see, for 
instance, [||]): The elementary charge excitations of the Hubbard model at finite 
positive U are spinless, hence the density is smaller by a factor of two. 

To complete our understanding of the gas phase of the one-dimensional Hub- 
bard model let us consider the low temperature thermodynamics of the gas phase. 
The thermodynamics of the Hubbard model was first considered by Takahashi 
[[Rp. The limiting case we are interested in, however, was not studied in Taka- 
hashi's paper. Takahashi expressed the Gibbs free energy co = — P (P pressure) in 
terms of the dressed energies K(k), £ n (A), £^(A), of elementary excitations at fi- 
nite temperature. k(&) is the dressed energy of particle (or hole) excitations, £„(A) 
describes spin excitations and £',(A.) so-called k-A strings [p3|]. All k-A strings are 
gapped [ |T4| . They do not contribute to the low-temperature thermodynamic prop- 
erties of the Hubbard model [ p~5] ] and drop out of the equation for the pressure, 
which simplifies to 

Similarly, the integral equations for the dressed energies at low temperature be- 
come 

oo 

k(£)=-,u-2cos£- T £ ([n]ln(l+e-TU (sinfc) , (12) 

n=l 

--— J dkcoska n (A—smk)\ny\+e T J 

oo 

+ £ (A, im ln(l+ e -T^( A ) , (13) 



In (l+e t 



where n = 1,2,3,... in equation (|T3|), and 



6 



[n] and A nm are integral operators defined by 

([0]/)(A)=/(A), (15) 

/CO 
dA'a n (A-A')f(A'),n=l,2,--- (16) 
-oo 

min{/?,m} 

A„ ra = £ ([|n-m|+2(7-l)] + [|n-m|+27]). (17) 

The gas phase is characterized by the absence of a Fermi surface for K(k) . Thus 
K(k) is positive in the zero temperature limit, and the first term on the right hand 
side of (JT3J) becomes exponentially small in T. Dropping this term, the equations 
( |13l ) decouple from (JT2|) . Since the equations become independent of A, it is not 
hard to solve them. The solution, exp{£ ;i (A)/r} = n{n + 2), is the same as in the 
infinite coupling limit U — > °° (cf. e.g. Inserting this solution into dl2| ) we 

obtain 



K(Jfc) = -A/-2cosfc-rin2 . (18) 

Our initial assumption, that limr^o K (^) > holds for all k, is self-consistent, if 
fj + 2 < 0, which is precisely the condition for being in the gas phase stated above. 



With fll8| ) the low temperature expression for the pressure becomes 



P=— dkln(l+ 2e t ) w W — e~ , (19) 



2% J -n V / V 7t 

and we see that the density D = dP/djji and the pressure P are related by the ideal 
gas law, 

P = TD . (20) 

There are two important lessons to learn from our simple calculation. First, the 
low temperature limit in the gas phase works the same way as the strong coupling 
limit at finite temperatures. Second, the low temperature Gibbs free energy CO = 
— P in the gas phase shows no signature of the discreteness of the lattice. It is 
the same as for the impenetrable electron gas (see below), which is a continuum 
model. This agrees well with our intuitive understanding of the gas phase at low 
temperature: (i) The mean free path (= \/D) of the electrons is large compared to 
the lattice spacing (which we set equal to unity so far), (ii) Their kinetic energy is 
of the order T . Hence, the effective repulsion is large for T <C t7. (iii) The ideal 
gas law holds at low temperature. 



7 



3 Scaling 



The above arguments show that only electrons with small momenta, correspond- 
ing to long wavelengths contribute to the low-temperature properties of the Hub- 
bard model in the gas phase. Thus the Hubbard model in the gas phase at low tem- 
perature is effectively described by its continuum limit. In order to perform the 
continuum limit we have to introduce the lattice spacing A and coordinates x = An 
connected with the nth lattice site. The total length of the system is £ = AL. The 
continuum limit is the limit A — * for fixed i. In this limit we obtain canonical 
field operators ^ Q {x) for electrons of spin o as 



T H = A 2 T , m + 2 = A 2 ij ,k H = Ak,t H = t/A 2 , B H = A 2 B , (22) 



where k denotes the momentum, t the time and B the magnetic field, which we 
shall incorporate below. The index refers to the Hubbard model. Then, in the 
limit A — > 0, we find 



Here H is the Hamiltonian for continuous electrons with delta interaction, 



r l l z r U 
H = j_ e/2 dx{(d x V+(x))d x V a (x) + -¥+(*)¥+(*)^(x)¥ T (*) 

-rt(jc)¥a(*)} • (24) 



Note that the coupling c\ = U / A of the continuum model goes to infinity! This 
is a peculiarity of the one-dimensional system. The effective interaction in the 
low density phase becomes large. Similar scaling arguments lead to an effective 
coupling C2 = U in two dimensions and to C3 = AU in three dimensions, i.e. unlike 
one-dimensional electrons three-dimensional electrons in the gas phase are free. 

4 Universality 

What happens to more general Hamiltonians in the continuum limit? Let us con- 
sider Hamiltonians of the form Hq = Hh + V, where Hh is the Hubbard Hamilto- 
nian and V contains additional short range interactions. We shall assume that V is 
a sum of local terms Vj which preserve the particle number. Then Vj contains as 




(21) 



Let us perform the rescaling 



H H /T H =H/T. 



(23) 
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many creation as annihilation operators, and the number of field operators in V; is 
even. We shall further assume that Vj is hermitian and space parity invariant. 

According to equation ( |2"T| ) every field cj, a on the lattice contributes a factor 
of A 1 / 2 in the continuum limit. One factor of A is absorbed by the volume ele- 
ment dx = A, when turning from summation to integration. Thus, if Vj contains 
8 or more field operators, then V ~ A 3 and V/Th vanishes. If Vj contains 6 field 
operators, then at least two of the creation operators and two of the annihilation 
operators must belong to different lattice sites, since otherwise Vj = 0. A typical 
term is, for instance, Vj = c ^t c n c 1+i | c 7+i.T c /,l c /,T- m me continuum limit we 
have c J+L j = A 1 / 2 ^ (jc) + A 3 / 2 3 X V I / | (x) + <9(A 5 / 2 ). Hence, the leading term van- 
ishes due to the Pauli principle. The next to leading term acquires an additional 
power of A. We conclude that V ~ A 3 and thus V/Th — > 0. 

If Vj contains 4 fields, then 

V ~ A 2x ¥^{x) x ¥+(x) x ¥ i {x) x ¥ 1 {x) + 0(A 4 ) . (25) 

Here the first term on the right hand side is the density-density interaction of 
the electron gas. In order to arrive at the impenetrable electron gas model the 
coefficient in front of this term has to be positive. Note that there are no terms of 
the order of A 3 on the right hand side of ( p5| ) and thus no other terms than the first 
one in the continuum limit. Terms of the order of A 3 would contain precisely one 
spatial derivative. They are ruled out, since they would break space parity. 

Considering the case, when Vj contains 2 fields, we find, except for the kinetic 
energy and the chemical potential term, terms which correspond to a coupling to 
an external magnetic field Bh- For these terms to be finite in the continuum limit 
we have to rescale the magnetic field as Bh — A 2 B (cf. equation (K2I)). 



Our considerations show that the impenetrable electron gas model with mag- 
netic field, 



ft/2 

H B =H+B J ^dx Wo^p (x) , (26) 



is indeed the universal model (for small T) for the gas phase of one-dimensional 
lattice electrons with repulsive short-range interaction. 



5 Impenetrable electrons 

The impenetrable electron gas is the infinite coupling limit of the electron gas with 
repulsive delta interaction (A — > in (^)), which was the first model solved by 
nested Bethe ansatz [jl], §]. The pressure of the system as a function of T, n and B 
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is known explicitly [p^|], 



*ln 1+ — +e— , (27) 



27C 

and may serve as thermodynamic potential. The expression ( |2"7| ) is formally the 
same as for a gas of free spinless Fermions with effective (temperature dependent) 
chemical potential /j e ff = /u + T ln(2 coshB/T) . Hence the Fermi surface vanishes 
for limr^o/Jeff = fx+\B\ < 0. The finite temperature correlation functions of the 
impenetrable electron gas depend crucially on the sign of fi e ff- This allows us 
to define the gas phase at finite temperature by the condition /j e ff < 0, which is 
also sufficient for deriving the ideal gas law ( |20"1 ) from the low temperature limit 
of (p7|). Note that for zero magnetic field and small temperature equation ( |27j ) 



coincides with the right hand side of (|I9|). 

The time and temperature dependent (two-point) Green functions are defined 

as 

tr^-^/ r * T (x,f)^f(0,0)' 
G U x i f ) = tr(e- H B/T) <28> 

trfe-fi^Wtfat)^ (0,0)' 

%M = ~ — ^mjf] -• (29) 

For the impenetrable electron gas these correlation functions were represented as 
determinants of Fredholm integral operators in JT7| , [T8| ]. The determinant repre- 
sentation provides a powerful tool to study their properties analytically. 

In [Q, ^ the determinant representation was used to derive a nonlinear par- 
tial differential equation for two classical auxiliary fields, which determine the 
correlation functions. This partial differential equation is closely related to the 
Heisenberg equation of the quantum Hamiltonian ([24]). It is called the separated 
nonlinear Schrodinger equation. Together with a corresponding Riemann-Hilbert 
problem it determines the large-time, long-distance asymptotics of the correla- 
tors (PS|), ( |2"9"| ) (for details see the following sections). In [Q, |5p the asymptotics 
x,t —> °° was calculated for fixed ratio £o =x/2t. The crucial parameter for the 
asymptotics is the average number of particles xD in the interval [0,x]. If x is 
large but xD 1 (i.e. T small), an electron propagates freely from to x, and the 



correlation functions <J28D, fl29p are those of free Fermions, 



G}{x,t) = 'rteWe*, (30) 



in t 

p ~~K {u-B-ki) . . . , :„2 

G - M = iLe^T^t-le-tMe-* . (31) 
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The true asymptotic region is characterized by a large number xD of particles in 

the interval [0,jc], specifically, xD > z~K where z c = (T 3 / 4 e- k o/ 2T )/(2n l / 4 k 3 Q /2 ). 
If the latter condition is satisfied, the correlation functions decay due to multiple 
scattering. The cases B > and B < have to be treated separately. For B > a 
critical line, x = 4t\^B, separates the x-t plane into a time and a space like regime. 
The asymptotics (for small T) in these respective regimes are: 

Time like regime (x < At\J~B): 



t THzc) e -xD [ 

^/AnzcxD i 



G^. (x, t) = Gf (x, t) 7==^== ' , (32) 



where 

v(Zc) 



2D,^ / V*o/2r ! fr (U+B]/T 
ft; , D\ = -\ -e^ +B)/T . (33) 

Dj = dP/d(/j + B) is the low temperature expression for the density of down-spin 
electrons. 

Space like regime (x > AtyB): 

G±(x,t) = Gf(x,t)t* iv ^e- xD l , (34) 



where 



e (3B+v-k 2 )/T 



> = s— ■ (35) 

For B < there is no distinction between time and space like regimes. The asymp- 
totics is given by (|3~4"1). 

It is fair to mention here that the calculation of the asymptotics (j3~2"|), (p4"|) is 



rather lengthy. Equations ( |32D and ( p4h are asymptotic expansions in f consisting 
of an exponential factor, a power law factor and a constant factor. Note that the 
method employed in allows for a systematic calculation of the next, subleading 
orders. 

The leading exponential factor in ( |32| ) and ([34]), has a clear physical inter- 
pretation: Because of the specific form of the infinite repulsion in (|Z4"|), up-spin 
electrons are only scattered by down-spin electrons. This is reflected in the fact 
that the correlation length is 1/Di. The expression l/D^ may be interpreted as 
the mean free path of the up-spin electrons. Thus the correlation length for up- 
spin electrons is equal to their mean free path. The exponential decay of the 
two-particle Green functions means that, due to the strong interaction, an up-spin 
electron is confined by the cloud of surrounding down-spin electrons. Thus we 
are facing an interesting situation: Although at small distances the electrons look 
like free Fermions, they are confined on a macroscopic scale set by the mean free 
path 1/D|. 
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6 Outline of the derivation 

The derivation of the above results on the asymptotics of the two-particle Green 
functions at low temperature is based on the fact that the impenetrable electron 
gas model is exactly solvable by Bethe ansatz. The Bethe ansatz eigenfunctions 
[jTJ, 0] and the thermodynamics of the model are known since long. But only 
recently a determinant representation for the two-particle Green functions was 
derived by Izergin and Pronko tf\J\, [TSp. Their derivation includes the following 
steps: 

(i) A change of basis for the spin part of the Bethe ansatz wave function from 
inhomogeneous XXX to XX spin chain eigenfunctions, which is possible at 
infinite repulsion. 

(ii) The calculation of form factors in the finite volume. 

(iii) A summation of the form factors. 

(iv) The thermodynamic limit. 



The details of the calculation can be found in the article Q18[]. 

The asymptotic analysis of the correlation functions was performed in [Q, Bp. 
Starting point was the determinant representation of Izergin and Pronko (sec- 
tion 0), which is valid for all x and t. The non-trivial ingredients of the deter- 
minant representation are certain auxiliary functions b ++ and B and the Fred- 
holm determinant det(/+ V ) of an integral operator V. A direct, yet lengthy cal- 
culation shows that b ++ and B satisfy the separated non-linear Schrodinger 

equation (section |]), which is a well-known integrable partial differential equa- 
tion. The logarithm of the Fredholm determinant plays the role of its tau-function 

(section |9|). Moreover, a Riemann-Hilbert problem that fixes b+ + and B as 

solutions of the separated non-linear Schrodinger equation can be derived from 
the determinant representation (section [Tj]). The Riemann-Hilbert problem is the 
appropriate starting point for the asymptotic analysis of the correlation functions 
Gi (section |TT|). 

Luckily, the differential equation and the Riemann-Hilbert problem turn out 
to be of the same form as in case of the impenetrable (spinless) Bose gas [ |T9[ 



20Q . Therefore a theorem obtained in the asymptotic analysis of the impenetrable 



Bose gas pT|] could be applied to the impenetrable electrons as well. In contrast 



to the bosonic case, there is, however, an additional external integration in the 
determinant representation of the impenetrable electron gas. This integration can 
be carried out in the low temperature limit, by the method of steepest descent 
(section |l2| ). 
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7 Determinant representation 

Let us recall the determinant representation for the correlation functions Gr£ (x, t) , 
which was derived in [jT7|, [T8p. We shall basically follow the account of [|]]. Yet, 
it turns out to be useful for further calculations to rescale the variables and the 
correlation functions. The rescaling 

x r = -Vfx/2 , t r = Tt/2, (36) 
g± = G±/VT, (37) 
$ = Meff/T , h = B/T (38) 

removes the explicit temperature dependence from all expressions. Furthermore, 
it will allow us to make close contact with results which were obtained for the 
impenetrable Bose gas [ |T9| , |2"T| , |2T)| ]. The index V in d3"6] ) stands for 'rescaled'. 



For the sake of simplicity we shall suppress this index in the following sections. 
We shall come back to physical space and time variables only in the last section, 
where we consider the low temperature limit. 

The rescaled correlation functions g + and g~ in the rescaled variables can be 
expressed as [Q, |5]], 

^■'^ s ■ (39) 

g -ai(P-A-ln(2ch(A))) f n 

g ~ M= 4^ J_dnF(y^)B—det(t+V) . 

Here yand F(y,r|) are elementary functions, 



(40) 



y = l+e 2/ \ (41) 

gill 

F(y,ti)=H — H ^. (42) 

det(7 + V) is the Fredholm determinant of the integral operator / + V, where / is 
the identity operator, and V is defined by its kernel V(k,fi). X and /u are complex 
variables, and the path of integration is the real axis. In order to define V(k,/u) we 
have to introduce certain auxiliary functions. Let us define 

x(X) =i(X 2 t + Xx) , (43) 

<W =TT^p ' (44) 

poo e -2x(fi) 

E(X)=p.v. d/u- ^, (45) 
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e-M=\l^eW t (46) 



e + {X) =^^e- x W{(l-cos(Ti))e 2T W J E(X)+sin(Ti)} . (47) 

Note that ~&(X) is the Fermi weight. V(X,/u) can be expressed in terms of e + 
and e_, 

y (M= ^(^-(A*)-^-(>.) . (4g) 

A — 

Denote the resolvent of V by 

(7 + V) (/-/?) = (7-7?) (7 + V) =/. (49) 
Then R is an integral operator with symmetric kernel [|19|], 

^ )= MMiMM), ,50) 

A, — /j 

which is of the same form as V(k,/j). The functions f± are obtained as the solu- 
tions of the integral equations 

/oo 
d l iV{\,n)f±{ i i) = e±{X). (51) 
-oo 

We may now define the potentials 

/oo poo 
dXe a (X)f b (X) , C ab = dXXe a (X)f b (X) (52) 
-oo J — oo 

for a,b = ±. B enters the definition of g~(x,t), equation (|40|). & ++ in ( |39| ) is 

defined as 

fc ++ =5-H--G(jc,0, (53) 

where 

The remaining potentials B ab and C ab will be needed later. 

It is instructive to compare the determinant representation ([39|) for the correla- 
tion function g + (xj) with the corresponding expression for impenetrable Bosons 
(cf e.g. page 345 of [|2"0|]). The main formal differences are the occurrence of 
the T] -integral in < p9] ) and the occurrence of T| in the definition of e+. As can be 
seen from the derivation of ( j3~9"t ) in [|18|], the T) -integration is related to the spin de- 
grees of freedom. Fort] = ±% the expression — je 2i $ t b++det(I + V) agrees with 
the field-field correlator for impenetrable Bosons (recall, however, the different 
physical meaning of (3). 
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8 Differential equations 

As in case of impenetrable Bosons JT9] , |ZD| ] it is possible to derive a set of inte- 

grable nonlinear partial differential equations for the potentials b++ and B and 

to express the logarithmic derivatives of the Fredholm determinant det(7 + V) in 
terms of the potentials B ab and C a t,. 

The functions f± satisfy linear differential equations with respect to the vari- 
ables x, t, and (3, 

L({ + )=M( f ;)=N( f ;)=0 : (55) 



J- J \f-J \f- 

The Lax operators L, M and N are given as 

L=d x + ilo z -2iQ 1 (56) 
M =d t + iX 2 cf - 2iXQ + d x U , (57) 

N =2X3 p + d x + 2itlo z + bco z - AitQ - 2d$U , (58) 
where the matrices Q and U are defined according to 

fl=r.° b ?) • VJ- B +- M. (59, 



5 J \~ B B 

Mutual compatibility of the linear differential equations ( |55| ) leads to a set 

of nonlinear partial differential equations for the potentials b ++ and B In 

particular, the space and time evolution is driven by the separated nonlinear Schro- 
dinger equation, 

id t b++ = - 4b 2 ++ B^ , (60) 

id t B_ = \d 2 x B—+AB 2 __b+ + . (61) 



9 Connection between Fredholm determinant and 
potentials 

To describe the correlation functions (E% and (ffDh one has to relate the Fredholm 



determinant det(/+ V) and the potentials B a t> and C a \,. Let us use the abbreviation 
a(x, t, (3) = lndet(7 +V) . The logarithmic derivatives of the Fredholm determinant 
with respect to x, t and [3 are 

d x c = - 2iB + , (62) 

dtC = -2i(C+-+C-+ + G(x,t)B—), (63) 

3pG = - 2i/3p(C+_ + C_ + + G(jc, t)B— ) - 2u3 p 5 + _ - 2(ap5 + _) 2 

-2if( J s__ap& ++ -^ ++ ap J s__)+2(a p Z7 ++ )(ap J B__) . (64) 
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For the calculation of the asymptotics of the Fredholm determinant we further 
need the second derivatives of o with respect to space and time, 

d 2 x o =AB__b ++ , (65) 
d x d t c =2i(B—dJ>++ - b++d x B—) , (66) 
dfo =2i(B—d t b++ - b++d t B—) + SB 2 __b 2 ++ + 2(d x B—)(d x ff + +) . (67) 

Note that 

lim o = . (68) 

This follows from limp^_ 00 f>(X) = and is important for fixing the integration 
constant in the calculation of the asymptotics of the determinant. 

10 The Riemann-Hilbert problem 

From now on we will restrict ourselves to the case of negative effective chemical 
potential, [3 < 0. Recall that this is the condition for the system to be in the gas 
phase. For negative (3 the logarithmic derivatives d x o and d t o of the Fredholm de- 
terminant and the potentials b ++ and B are determined by the following matrix 

Riemann-Hilbert problem, which was derived from the determinant representation 
(see section 0) in @] . 

(i) 4> : C -> End(C 2 ) is analytic in C \ R. 

(ii) limj^^X) =I 2 . 

(iii) (|) has a discontinuity across the real axis described by the condition 

*-w=*+wG w ^> (fmZi) (69) 

for all XgM. 

Here I 2 denotes the 2x2 unit matrix. The functions p(X) and q(k) are defined as 
p(X) =i(cos(r|) - 1)(1 -$(\))a+(X)a-(\) , (70) 

a + (A)a_(A) 



where 



a w=-»{^/-"^ ln(1+ * ( " )(e "- 1)) } ' (72) 
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The functions d x G, d t G, b ++ and B can be expressed through the coeffi- 
cients in the asymptotic expansions of §(X) and ln(a(X)) for large spectral pa- 
rameter X. Let 

(b« d»( 2 ) / 1 \ 

$(%)=I 2 + V— + \- + „ (73) 



X X 2 \ V 

and 



ta («W) = T + i + (pl- < 74 > 



Then 



a x a = 2iai+itr{4>( 1 )o z } , = 4ia 2 + 2itr{4>( 2 )a z } , (75) 
£ ++ =(^ , 5__ = -4 1 1 ) . (76) 

The Riemann-Hilbert problem is the appropriate starting point for the asymp- 
totic analysis of the potentials b ++ and B which determine the asymptotics of 

the two-particle Green functions G*+. For impenetrable Bosons a similar analy- 
sis was carried out in pT|]. Fortunately, the result of [ |2~T| ] depends only on some 
general properties of the functions p(X) and q(X) entering the conjugation matrix 
in (j69|), and also applies in the present case. Alternatively, the non-linear steepest 



descent method of Deift and Zhou p2j] could be applied. 



11 Asymptotics of the correlation functions 

The direct asymptotic analysis of the Riemann-Hilbert problem yields the leading 
order asymptotics (x, t — ► °° for fixed ratio Xq = —2x/t) of the functions d x o, d t o, 

b ++ and B [|5|, |2l|]. It turns out, in particular, that b ++ and B are a decaying 

solution of the separated nonlinear Schrodinger equation (PH), (|6"T|). Now the 
form of the complete asymptotic decomposition of the decaying solutions of the 
separated nonlinear Schrodinger equation is known p3] , (Z4[ ] . 

_i / - 2» \n h At \ (be 2 . , „ 1 
6 ++ =? 2 I Mo +^ ^ __ M;;yt 1 ex p|— -ivln4?| , (77) 

i / 00 2 " \n k At \ ( \x 2 ^\ 
B U + £ £ v,* exp +ivln4r , (78) 

where wo, vo, w n &, and v are functions of Xq = —x/2t and of (3 and r). Inserting 
the asymptotic expansions for B and b ++ into the differential equations (|5H|), 
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( PI ) we obtain expressions for u n k, v n k and v in terms of uq and vq, i.e. the two un- 
known functions uq and vo determine the whole asymptotic expansion (f77|), (]78|). 
But mo and vo are obtained from the asymptotic analysis of the Riemann-Hilbert 
problem (for the explicit expressions see [Q]). Hence we know, in principle, the 

complete asymptotic decomposition of the potentials b ++ and B 

In order to obtain the asymptotics of the two-particle Green functions we still 
need the asymptotics of the Fredholm determinant. The Fredholm determinant 
is related to b ++ and B through equations (|j5|)-(]57[) and (p3)-(p^). We may 



integrate (|65t)-(|67l) to obtain the asymptotic expansions of d x o and d t o. The inte- 
gration constant is a function of (3. It is fixed by the leading asymptotics, which, 
using (f75|), can be obtained from the direct asymptotic analysis of the Riemann- 
Hilbert problem. Then, integrating (|o"2|)-(|53|) yields a up to a numerical constant, 
which follows from the asymptotic condition (|68|). The calculation is the same as 
for the impenetrable Bose gas and can be found on pages 455-457 of [E0]. 

Finally, we obtain the following expressions for the leading asymptotics of the 
correlation function, 



+ (x,t) =e a / 2t + 2it ^ e - 2it ( h +H^h(h))) I dr ^ 



71 



F(y,t\) 



l-Ti l-cos(ri) 

•C+(X ,p,ri)(40^ v - l)2 exp|iy' OO ^|x + 2^|ln((p(?i,p))| , (79) 



(jC,t) =e -« 2 A-2kP e 2 1 r(/ 1 +ln(2ch(/ ; ))) r dr[ 

J— % 



2y 



• C (Xo, p, ti) (40^ (v+i)2 exp ( - r dl \x + 2lt\ ln((p(A, (3)) 

I Tl J-oo 



(80) 



where 



<p(X,P) =1 ^- iT lsign(X-Xo) _ X \ (81) 

v = -^-ln(l-2(l-cos(TiMXo)(l-^o))) , (82) 
In 

C + (Ao, P, tl) = - I Binft/2) | ^ exp 1 1 (A* - p) + i^o + y 

-/ Jp(iv/2 + va p ^ ) (83) 
+^2 f_J$ ( 3 P /^sign(A-A )ln((p(A,p)) 

C-(Ao, P,ri) =C+(A , p,ri) exp(-(Aj - p) - 2i¥ )/ sin 2 (ri/2) . (84) 
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A,o = —x/2t is the stationary point of the phase x(k) (i.e. x'(A,o) = 0), and the 
functions and *Fi are defined as 

371 

'P =- T +aigr(iv)+'Pi J (85) 
JXsign(X-Xo)ln|X-X |3^1n((p(X,(3)) . (86) 

Equations ( J79| ) and ( |S0| ) are valid for large f and fixed finite ratio Xq = —x/2t. 



Correlations in the pure space direction t = were discussed in Q25j]. We would 
like to emphasize that ([79]) and (|8U|) still hold for arbitrary temperatures. The low 
temperature limit will be discussed in the next section. Note that there is no pole 
of the integrand at r\ = 0, since a/v ~ |r|| for small r\ and thus C + (Xq, p\r|) ~ r\ 2 . 



12 Asymptotics in the low temperature limit 

For the following steepest descent calculation we transform the T]-integrals in (f79|), 
( PU[ ) into complex contour integrals over the the unit circle, setting z = e" 1 . Since 
we would like to consider low temperatures, we have to restore the explicit tem- 
perature dependence by scaling back to the physical space and time variables x and 
t and to the physical correlation functions G^j. Recall that in the previous sections 
we have suppressed an index V referring to 'rescaled'. Let us restore this index 
in order to define k = l Vf = x/2t, -&(k) = -& r (k/Vf), (p(&, (3) = <p r (k/y/T, (3), 
C*(*o,P,z) = C*(Ao,M),F(y,z) =F r (y,Ti). Then 

F(Y,z) 



Mv(z)-iy e ts( z ) 



(87) 

G^(x,t) =-iVfe-' vc2 / 4t - [t ^ jdz^^-C-{kQ^z){2Tt)± 



(v(z)+i) 2 e tS(z) 



(88) 



where 



S( z ) = - f dk\k-k \ In(q>(Jfe, (3)) . (89) 

71 J-oo 



We would like to calculate the contour integrals (|87|), (]88|) by the method of 
steepest descent. For this purpose we have to consider the analytic properties of 
the integrands. Let us assume that ko > 0, and let us cut the complex plane along 
the real axis from — °° to — and from — gP~ fc o/ r to 0. The integrands in ( |8"7j ) 
and <J8"%D can be analytically continued as functions of z into the cut plane with the 
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only exception of the two simple poles of F(y,z) at z = f . We may therefore 
deform the contour of integration as long as we never cross the cuts and take into 
account the pole contributions, if we cross z = J or z = Y _1 . 

The saddle point equation dS/dz = can be represented in the form 

die k (' die Jc 

(90) 



i +z -i e -$ e (k-k y-/T j i +ze -^ e (k+k l) )yT ■ 

This equation was discussed in the appendix of [Q]. In [|]] it was shown that ( |90D 
has exactly one real positive solution which is located in the interval [0, 1]. It was 
argued that this solution gives the leading saddle point contribution to ( |S7| ) and 



(|88|). At small temperatures (gOJ) can be solved explicitly. There are two solutions 
z± = ±z c , where 

2% l / 4 kl /2 

In the derivation of ( ]9T1 ) we assumed that ko 7^ 0- The case ko = has to be treated 
separately (see below). 

The phase tS(z) has the low temperature approximation 

tS(z) = -2k Dt I M - i J i c + 1 - z\ . (92) 

Here D = dP/d/jis the density of the electron gas. The low temperature expansion 
( p2| ) is valid in an annulus e$~ k °/ T \z\ e~$, which lies in our cut plane. The 
unit circle and the circle \z\ = z c are inside this annulus. We may thus first apply 
( p2"| ) and then deform the contour of integration from the unit circle to the small 
circle \z\ = z c - Let us parameterize the small circle as z = z c e m , oc E [—71, 7c] . Then 
S(z(ct)) = — 2koD((z c — l) 2 + 2z c (l — cos(a))), which implies that the small circle 
is a steepest descent contour and that on this contour S(z-) < S(z) < S(z+). The 
maximum of S(z) on the steepest descent contour at z = z+ is unique and therefore 
provides the leading saddle point contribution to (]87|), ( |88| ) as t — > oo. The saddle 
point approximation becomes good when tS(z(o£)) = —2koDt((z c — l) 2 +z c a 2 + 
O (oc 4 )) becomes sharply peaked around a = 0. Hence, the relevant parameter for 
the calculation of the asymptotics of G** is Ik^Dt = xD rather than t. xD has to 

be large compared to z^ 1 ■ The parameter xD has a simple interpretation. It is the 
average number of particles in the interval [0,*]. Let us consider two different 
limiting cases. 

(i) xD — > 0, the number of electrons in the interval [0,x] vanishes. In this regime 
the interaction of the electrons is negligible. An electron propagates freely 
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from to x. G TT cannot be calculated by the method of steepest descent. 



TT 



We have to use the integral representation (f79|), ( J80| ) instead. Since tS(z) 
and v(z) tend to zero on the contour of integration, the integrals in < J79| ) and 

(see ©, O), which is the 



are easily calculated. We find G 
well known result for free Fermions 



± 
11 



(ii) xD ^> Zc 1 , the average number of electrons in the interval [0, x] is large. This 
is the true asymptotic region, x — > °°. In this region the interaction becomes 
important. At the same time the method of steepest descent can be used to 
calculate G^j . This case will be studied below. 

In the process of deformation of the contour from the unit circle to the small 
circle of radius z c we may cross the pole of the function F(y, z) at z = Y~ 1 . Then we 
obtain a contribution of the pole to the asymptotics of G^. It turns out that the pole 

contributes to G++, when the magnetic field is below a critical positive value, B c = 
&q/4. Below this value the pole contribution always dominates the contribution of 
the saddle point. Hence, we have to distinguish two different asymptotic regions, 
B > B c and B < B c . On the other hand, if we consider the asymptotics for fixed 
magnetic field, we have to treat the cases B > and B < separately. For B > we 
have to distinguish between a time like regime (&g < 45) and a space like regime 



(&q > 45). In these respective regimes we obtain the asymptotics (p2|), (]34[)- 



In the limit B — > — °°, p — > — oo, p — B fixed there are no |-spin electrons left 
in the system, Z)j — > 0. This is the free Fermion limit. In the free Fermion limit 
B < B c , and the asymptotics of GfJxj) and G7Jx,t) are given by the equations 



([34[), which turn into the expressions (p0|), ( [31] ) for free Fermions. 

The pure time direction ko = requires a separate calculation. For ko = the 
saddle point equation ( ]9(1| ) has the solutions z = ±1 for all temperatures. The unit 
circle is a steepest descent contour with unique maximum of S(z) atz— 1, which 
gives the leading asymptotic contribution to the integrals in (|7]) and (|SS|). We find 
algebraically decaying correlations, 



G|j (0, 



C + f 



G|! (0, ?) 



C Q t 



-l-itifi-B) 



(93) 



where 



:(l+2e- 2fi / r ) 



2V2nf 

( e (»+B)/T + e (M-B)/T^ ( { + e (v+B)/T + ^-B)/^ 



e 1 * 



l+2e- 2B / T 

'iV^T l+e 2B / T 



e (p+B)/T + e (p-B)/T 
l +e brt-B)/T +e 0i-B)/T 



(94) 



(95) 



These formulae are valid at any temperature. 
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